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In human societies opinion formation is mediated by social interactions, consequently taking place 
on a network of relationships and at the same time influencing the structure of the network and its 
evolution. To investigate this coevolution of opinions and social interaction structure we develop a 
dynamic agent-based network model, by taking into account short range interactions like discussions 
between individuals, long range interactions like a sense for overall mood modulated by the attitudes 
of individuals, and external field corresponding to outside influence. Moreover, individual biases can 
be naturally taken into account. In addition the model includes the opinion dependent link-rewiring 
scheme to describe network topology coevolution with a slower time scale than that of the opinion 
formation. With this model comprehensive numerical simulations and mean field calculations have 
been carried out and they show the importance of the separation between fast and slow time scales 
resulting in the network to organize as well-connected small communities of agents with the same 
opinion. 

PACS numbers: 64.60.aq, 87.23.Ge, 89.75.Fb 



I. INTRODUCTION 

The network approach has contributed significantly to 
our understanding of the structure, function, and re- 
sponse of various complex systems, from genetic tran- 
scriptions to human societies [3, [3] • In the case of human 
social systems, this approach was introduced by social 
scientists and they established important concepts and 
tools to study them [3, 0| • According to the social net- 
work paradigm social life consists of the flow and ex- 
change of norms, values, ideas, and other social and cul- 
tural resources channeled through a network [j|. Hence 
the networks serve as substrates for various collective 
social phenomena such as diffusion and spreading pro- 
cesses (of e.g. news and epidemics), opinion formation, 
language evolution, etc. 

Until recently the studied empirical data sets of so- 
cial systems remained rather limited since the basic 
source of data were questionnaires, thus the focus had 
been on small scale properties. However, the most re- 
cent development in information-communication technol- 
ogy has opened the possibility to collect much larger 
data sets from Internet, emails, phone records, etc. 
B 0, 0, H El E3- While the scope of these records 
is narrow as compared to detailed questionnaires, their 
huge amount and objective quantifiability enable us to 
study problems, which were impossible to treat before, 
including the investigation of the structure and dynamics 
of entire populations [l2| . Much has been learned from 
these studies such as the broad distributions of network 
characteristics, the small world properties, the modular 
organization of the social network, and the relationship 
between the network topology and the intensity of the 
ties in the net. 



One of the key problems in studying society and social 
dynamics falls under the term of "opinion formation" , 
which deals with a (measurable) response of the society 
to an issue, such as an answer to a political question or 
the acceptance of innovation. This problem has attracted 
a lot of interest and a number of models of competing op- 
tions have been introduced to study it [HI . The simplest 
one is the so-called voter model which has a binary 
opinion variable with the opinion alignment proceeding 
by a random choice of neighbors. Other discrete opin- 
ion formation models include those by Weidlich [l5| and 
Sznajd [16j . where more than just a pair of spins is associ- 
ated with the decision making procedure. These models 
are reminiscent to the Ising model and in the noiseless 
(T = 0) case they usually approach an asymptotic state 
of consensus (ferromagnetic state). 

If these models are studied on complex networks in- 
stead of regular lattices [13, EH, Ell, the temporal dy- 
namics of the phase transition from one type of con- 
sensus to the other one may change significantly. In- 
troducing noise or "social temperature" makes the anal- 
ogy with physics even closer and leads to a number of 
new phenomena from paramagnetic phase to stochastic 
resonance. In addition the introduction of an external 
field as the carrier of the influence of mass media is a 
natural and widely considered generalization (see, e.g., 
[13, H3])- Also systems with more than two possible opin- 
ions have been considered [2l| , which have found natural 
extension in models where the opinion of individuals is 
represented by a continuous variable [13, H3, [24|. By 
introducing a parameter of bounded confidence or tol- 
erance that describes the situation where opinions close 
enough approach each other, the asymptotic state can 
be quite rich in allowing the coexistence of a number of 
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different opinions. We also note that opinion dynamics 
is closely related to problems of competing cultures (2f| 
or languages [261 ], 

Opinion formation in a human society is mediated by 
social interactions between individuals, consequently it 
takes place on a network of relationships and at the same 
time it influences the network itself. Such behavior was 
analyzed in a seminal paper by Holme and Newman [27| . 
Recently this problem of coevolving opinion formation 
networks was considered for discrete [IM and contin- 
uous [3(3] opinion dynamics. In these the separation of 
the time scales was controlled by the relative frequency of 
the opinion updating and rewiring processes, which is an 
important parameter of the problem. It was found that 
the adjustment of the network has a major effect both on 
the opinion dynamics and on the resulting topology and 
that minor changes in the updating rules have dramatic 
consequences. In the case of discrete models it turned 
out that the self-adaptation of the network resulted in 
the symmetry breaking of the asymptotic state with re- 
spect to the direction of the information transfer. In the 
case of continuous opinion formation dynamics rewiring 
tends on one hand to hinder the consensus formation for 
large tolerance since it breaks the clusters. On the other 
hand there is a tendency to decrease the number of differ- 
ent clusters for small tolerance. In all these investigations 
agents were considered to be uniform. 

Here we also consider the continuous opinion forma- 
tion in a coevolving network but with the addition that 
the individuals need not to be uniform. In particular we 
examine a situation in which these non-identical individ- 
uals form their opinions in information-transferring inter- 
actions with others. This could happen directly through 
discussions between individuals constituting the network 
structure or indirectly by sensing the overall mood or 
opinion of all the other individuals, the effect of which 
would depend on their personal attitude towards the 
overall opinion. We expect that in general the time scale 
for detectable changes in the network structure is con- 
siderably slower than the time scale for the direct social 
exchanges to take place. In order to describe this situa- 
tion we have developed a dynamic network model, where 
we take into account short range interactions for direct 
discussions between pairs of individuals, long range inter- 
actions for sensing the overall opinion modulated by the 
attitude of an individual, and external field for outside 
influence. The opinion formation dynamics is assumed to 
take place with a fast time scale well separated from the 
slow time scale network topology coevolution included 
by opinion dependent link-rewiring processes of the kind 
introduced as the basic mechanisms of network sociology 
[3ll . HH . Note that unlike the model introduced in [3(| 
our model does not contain any parameter for bounded 
confidence, which could alter the time scale of network 
coevolution, as we discuss later. 

This paper is organized such that next we describe 
our model, followed by presenting the numerical results 
and then by finding analytical solutions to the dynami- 



cal equation with mean field approximation. Finally we 
discuss the results and draw conclusions. 



II. DESCRIPTION OF THE MODEL 

Let us consider an opinion formation in a network of 
a fixed number of individuals or agents (N) to whom a 
simple question is posed. We assume that the agents are 
unbiased but they have initial opinions concerning the 
issue of the question. The process of opinion formation 
is started by letting the agents to discuss or exchange 
their views with their acquaintances or friends for some 
time. After this some discussions are interrupted at reg- 
ular intervals to let discussions with new agents to start. 
At any time an agent may decide to fix its opinion to 
either total agreement or disagreement, after which the 
decision is considered irrevocable and no longer modifi- 
able, although the agent continues to discuss and influ- 
ence others. The whole process stops when everyone has 
balloted. It should be noted that here the aim is not to 
reach a consensus of opinion, but to give every agent a 
chance to use its social network to form a point of view 
through exchange of information. 

This kind of a social system can be described by a net- 
work where the nodes correspond to the agents and the 
edges to the links or interactions between them. Each 
individual agent i is described with the time dependent 
state variable x; £ [—1, 1] measuring the agent's instan- 
taneous inclination towards the posed question. The 
agent's interactions with other agents are described with 
the adjacency matrix elements Aij £ {1,0} that repre- 
sent the presence or absence of discussions between the 
pair of agents related to the question at hand. Thus the 
initial topology of the network defined by the adjacency 
matrix A changes when the agents are allowed to stop 
discussions. 

As for the evolution dynamics of this system we assume 
that the state variable Xi of node i depends on the links of 
the network, and likewise the links of the network change 
according to the values of the state variable of each node. 
These mutually dependent processes can be described in 
general with a dynamical equation of the following form 

^ = ^ +^d{x l ,x j ,g)A ij , (1) 
j 

which is based on the assumption of two separate time 
scales. Here we assume that many discussions take place 
before a change in the network structure occurs, that 
is, the time scale for each discussion or exchange of in- 
formation ("transaction") is dt, while the time scale for 
a change of connections in the network ( "generation" ) 
is T = gdt. The quantity g defines the number of 
transactions per generation and describes the separation 
between the fast transaction and slow generation time 
scales. In this equation O stands for an operator that 
changes the entries and/or the size of the adjacency ma- 
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trix. The actions of the operator O are necessarily dis- 
crete and limited to few processes. In general there are 
only four basic ways to modify a simple network: ei- 
ther one deletes/creates links or deletes/creates nodes. 
Then a common operation of rewiring can be regarded 
as a composed operator with one deletion followed by 
one creation of a link. An important point in modeling 
a particular system is that the operator O should con- 
tain rules for link deletion or creation reflecting all the 
additional information about the system. 

For large g or T the effect of discreteness of network 
evolution events on the agent's state variable Xi can be 
treated separately and the dynamics between such up- 
dates is essentially continuous. Then Eq. (fT]) can be ap- 
proximated as follows 

dx ' 

-q£ = otif {{xj}) + fi({xj})xi + h, (2) 

where the first term on the right hand side represents the 
combined effect of an agent i sensing the overall opinion 
of all the other agents (fo({ x j})) but modulated by the 
agent's own attitude towards overall or public opinion. 
We denote this personal attitude by on and assume it 
being random and uniformly distributed between -1 and 
1, where the former corresponds to completely oppos- 
ing and the latter to completely agreeing attitude of the 
agent i towards the overall opinion. The second term on 
the right hand side represents the direct discussions of the 
agent i with the agents j it is linked to. Note that the 
agents could be different also in their attitudes towards 
direct discussions, which would introduce another atti- 
tude parameter in the second term. Instead on should 
be considered as the relative attitude parameter. The 
third term on the right hand side, hi, is an external field 
representing the personal bias towards either opinion (-1 
or +1) due to e.g. mass media (newspapers, TV, radio). 
In Eq. ^ /o and /i standing for the dynamic long and 
short range interaction terms are defined as follows 

fo=Ej E *i. (3) 

and 

/l = sign(xj) x h ( 4 ) 

jemi(i) 

where mi{i) means the set of nodes that are t steps away 
from node i (or £ th neighbors of i) and £ m ax is the number 
of steps needed to reach the most distant neighbors of i 
(or the maximum range of interactions for i). 

Note that Eq. (|2|) allows changes of sign in the state 
variable but also its exponential decay or growth such 
that \xi\ could eventually become larger than one, which 
carries no meaning in our model. Thus we need first to 
detect the totally convinced agents and stop their dynam- 
ics, which can simply be done by a line of computer pro- 
gram "if abs(xi) > 1, then Xi — sign(xi)", that is, these 



agents cannot modify their state in subsequent times, but 
they are still linked to the network and taken into account 
in the dynamical evolution of the undecided agents. 

In order to consider the other dynamical process (and 
the corresponding time scale) contained in Eq. (fTJ), 
namely network connection topology changes at gener- 
ation time intervals T, we have here adopted the scheme 
of rewiring. There are basically two possibilities, either 
global or local rewiring, which we will perform with prob- 
ability y and 1 — y, respectively. These rewiring schemes 
are rooted to the fundamental link formation mechanisms 
of network sociology proposed by Kossinet and Watts 
(3l| : focal closure independently of distance thus being 
global and triadic closure between close network neigh- 
bors thus being local in nature. In both these schemes, 
first the agent i can choose to cut an existing link with an 
agent j, i.e. end a discussion if their opinions are incom- 
patible. In order to perform this process, the quantity 

Pij = Aij- - — (5) 

is calculated and all the links are put in ranking order. 
Then links with larger weights are deleted first since 
they correspond to divergence of opinion. After this link 
deletion step follows the link creation step using either 
local or global rewiring scheme. 

In the local rewiring scheme, an agent i can create a 
link with a second neighbor by starting a discussion with 
the "friend of a friend" if this new link can help the agent 
in reaching a state of total conviction (\xi\ = 1). In order 
to determine this we calculate the quantity 

gy = (l-^.)e[(A%] ^±^, (6) 

where the first two factors including the adjacency ma- 
trix (with Q[x] being the discrete step function Q[x] = 
J2k=i $k,x), test the existence of a link to the second 
neighbor and the third factor \xi + Xj\/2 is the measure 
of similarity of opinions between the agents. Then all the 
potential local rewiring links are put to ranking order, of 
which the links with larger weights qij are created first. 

In the global rewiring scheme, an agent i can create a 
link with further neighbors (£ > 2) provided their opin- 
ions are similar. In order to determine this we follow 
a similar procedure as above by calculating instead the 
following quantity 

r« = (1 - A^) (1 - 6 [(A% 3 \) (l JfLlfd) , (7 ) 

where once again the first two factors test the far link 
existence and the third the similarity of opinions between 
agents, followed by putting the potential links to ranking 
order and creating links with larger weights first. 

The rewiring of each node i is done in such a fashion 
that the number of links deleted is equal to the number of 
links created. Therefore, the total number of links should 
be conserved. This is clearly the case when the rewiring 
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is performed sequentially, but in case of parallel rewiring 
(as we do) it could happen that if both agents i and j are 
performing rewiring simultaneously, their mutual actions 
could lead to the net creation or deletion of a link. Then 
the total number of links would not be conserved. 

Next we will present the numerical results of our model 
followed by the analytical mean field theory treatment of 
the dynamical equation. 



III. NUMERICAL RESULTS 

For computer simulations we have first initialized our 
model system to a random network configuration of N 
nodes and average degree {k ). This is done in the begin- 
ning of each simulation run to secure a different random 
configuration for the initial network and good statistics 
for the averaged quantities of interest. As other initial 
conditions we chose the fixed constant Ui for the personal 
attitude and the state variable Xi(0) of agent i randomly 
from a uniform distribution between -1 and 1 and from 
a Gaussian distribution with unit standard deviation, re- 
spectively. In addition we chose all the external bias field 
terms to be hi = 0. 

With these initial choices we are left with only two 
parameters to vary, namely y and g, upon which the av- 
erage properties of the final network should depend and 
scale with the initial conditions (kg) and N. For the 
sake of simplicity we here set y = in all our numerical 
simulations, which means that we take into account only 
the local rewiring to study the effect of the parameter g 
alone. 

Since our model includes fast transaction dynamics 
and slower generation dynamics for network rewiring the 
simulations have been carried out by using the following 
two-step process. In the first step the dynamics of trans- 
actions described with Eq. (0) is realized by numerical 
integration using a simple Euler method, in which the 
time step was set to dt = 10 -4 as found to guarantee 
the stability and reliability of the numerical calculations. 
By keeping the parameters fixed the system is driven un- 
til the specified time (gdt) or g time steps to then do 
the second step, namely the network rewiring with the 
procedure described in the previous section. This two- 
step process is iterated until the system reaches its final 
state, where no more changes in the xi 's and in the net- 
work connections take place. We have found that the 
way the system approaches its final state consists of two 
regimes, which can be analyzed by controlling the num- 
ber of changed links per cycle. The first regime shows 
a rapid roughly exponential decay where both the num- 
ber of agents that change opinions and the number of 
rewired links is high. This fast decay crosses over to a 
very slow regime after most of the agents have reached 
either opinion x = 1 or — 1. During this slow decay, when 
frustrated links try to get optimal positions and the re- 
maining few undecided agents converge their opinions to 
x = 1 or — 1, the network structure is found not to change 



significantly. These few "non-conformists" that take ex- 
tremely long time to get fully convinced may serve as 
kind of nucleation centers if an external field (media ef- 
fect) is switched on, which is an issue we will investigate 
in the future. 

From the simulation results we calculated the following 
averaged single site properties: the degree (fc), the short- 
est path (L), the average clustering coefficient (C), the 
mean number of second neighbors (n 1 ^ ) , and the average 
cluster size or susceptibility (s) = n s s 2 / ^ s n s s that 
is the second moment of the number of s-size clusters, see 
e.g. [HI]. According to percolation theory the sums run 
over all the cluster sizes excluding the giant component, 
as identified in finite samples with the largest cluster. 

Fig.[T]shows the various properties of the final network 
as a function of g, averaged over 100 realizations for the 
initial average degree (ko) = 4 and for the system size 
N = 200. Here one can visually distinguish three regions 
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9 = T / dt 

FIG. 1: Average degree (k), average susceptibility (s), av- 
erage number of second neighbors ^n 1 - 2 ^, average shortest 

path (L) , and average clustering coefficient (C) for a network 
of size N — 200 and averages taken over 100 runs. The arrows 
indicate the g values where the susceptibility changes drasti- 
cally. The inset shows such quantities as a function of system 
size TV for g — 1000 and averaged over 1000 realizations. The 
graphs at the top are examples of final network configurations 
for N = 400 and g = 5, 10 3 , 10 s , which show three different 
types of configurations emerging. White circles mean agents 
with x — 1 and black those with x = — 1. 



of behaviour. For small g < 10 the dynamics of the 
state variable is practically irrelevant and accordingly, 



5 



the rewiring of the network is random and favours the 
formation of triangles. The final network appears divided 
into two large cliques of about N/2 nodes with almost 
no connections between them, as expected from models 
without dynamics. 

On the other hand, if g is larger than the mean number 
of iterations (i.e. > 10 4 ) needed for the state variable to 
reach either of its extreme values at all sites, no rewiring 
of the network is seen to occur. Thus the resulting net- 
work shows the properties of an Erdos-Renyi graph, i.e. 
(k) = N (C), and (L) = lniV/ln(fc}. Furthermore, since 
the initial network is random no apparent clustering of 
the two opposite opinions is seen to develop. 




FIG. 2: Communities found by the local fitness algorithm in 
the network configuration shown in the middle top of Fig. Q] 
for a system with N = 400, (fc ) = 4 and g = 1000. The com- 
munities are distinguished by arbitrary coloring, regardless of 
the value of variable x or opinion. 



For the intermediate g values both dynamical pro- 
cesses, slow and fast, are at play and the net result is 
that the final configuration of the network shows a segre- 
gation of several well-connected communities of roughly 
similar size and of uniform opinion, as the ones shown in 
the middle top of Fig. [TJ Such result makes this model 
of opinion formation suitable to describe the complex 
community structure often found in human societies (see, 



e.g. , [THHH). It is relevant to mention that Mandra et al. 
[35| have studied rewiring in the Ising-Glauber dynam- 
ics with tunable thresholds and identified the emergent 
domains of spins with communities. 

Such an emergence of medium size communities is 
clearly seen in the susceptibility (s), which peaks strongly 
(by two orders of magnitude) for the range of g from 
about 10 to more than 10 3 . It should be noted here 
that computationally the algorithm to solve Eq. (fJJ) in 
a rewiring network goes roughly as 0(N 3 ), which makes 
the simulation of systems with more than 1000 nodes 
very time-consuming. However, as it is seen in the inset 
of Fig. [TJ the system size dependence of all quantities is 
monotonic and can thus be extrapolated to the thermo- 
dynamic limit. Hence a subject of further study would be 
to investigate the nature of the transitions between the 
observed morphological phases. In addition, we would 
like to point out that all the other quantities show grad- 
ual behavior as a function of g. In the models of opinion 
formation without network dynamics, the resulting net- 
work is either a single cluster of consensus or two clusters 
of about the same size with opposite opinions. (As stated 
before we do not have bounded confidence in our model, 
which could be another source of the creation of many 
communities [29|.) 

Note that to produce the graphs in Fig. [T] we have 
used a visualization software [361 ] that tends to group 
agents according to their topological properties, and con- 
sequently communities are distinguishable by the eye. It 
is worth pointing out that within each community all its 
agents share the same opinion. In most previous studies 
the similarity of attributes between agents has not been 
considered to define communities (see for instance the 
recent comprehensive review by Fortunato [37|). 

In order to find out whether the communities formed 
by agents of the same opinion match the ones formed 
by topological similarities, we have used the community 
detection algorithm based on the local fitness concept de- 
scribed by Lancichinetti et al. [H| . In Fig. [5] we show 
the results for the network configuration at the middle 
top of Fig. Q] by using Eq. (1) of Ref. [H[, with the res- 
olution parameter set to 0.9. Observe that in Fig. Q] the 
color coding for the opinions allows the eye to distin- 
guish four well-defined communities, while in Fig. [2] the 
algorithm based on the local fitness separates the net- 
work into seven communities. Note that although there 
is in general good correspondence between these two par- 
titions, the community detection algorithm reveals sub- 
structure that the eye does not capture. 

The key point is to investigate what are the circum- 
stances that drive the system to form many communi- 
ties of approximately the same size. Here the role of 
the personal attitude parameter on is very important, 
since a network with agents having all positive a's should 
reach consensus quite rapidly without community struc- 
ture forming, while a network with only negative a's 
should separate into two clusters of different opinions. 
However, in case of a's being randomly distributed (as 
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is assumed in this study) frustration is introduced to the 
system, making the network structured. 

It is instructive to look at the distribution of the a's 
and how it relates to the cluster structure. For small 
values of g, where the rewiring process is very rapid and 
only two clusters eventually develop, the attitude param- 
eters have a minor role, so the a distributions in the two 
clusters are broad and similar. However, for the interme- 
diate g values the situation is different. Here the smaller 
clusters have a rather narrow distribution with mostly 
negative a values, while the distributions for larger clus- 
ters are broad and shifted toward positive a values. Nat- 
urally, the agents with negative a's do not feel comfort- 
able in a large homogenous cluster, thus they tend to 
build smaller communities. It is interesting that the per- 
sonal attitude of an agent has such an important effect on 
the clustering properties of the system in spite of the fact 
that it is not defined as a bias towards a certain opinion 
but rather as an agent's ability or intention to adjust or 
not to any kind of opinion of the other agents. 

The high connectivity and clustering coefficient of the 
final network configuration is the result of the appear- 
ance of communities: if there are X communities of size 
n, then N = Xn and if the clusters continue to be ran- 
dom, then (k) / (C) = n. Therefore, the region of values 
of g for which (k) / (C) is constant is where there are 
communities of size N/X. In Fig. [3] we show this ratio 
as a function of g for networks of different size. Notice 
that the region where communities are formed is wider 
for larger networks and that as the network size is in- 
creased the number of clusters X increases as well. In 
the inset we show the value of X for networks of var- 
ious sizes, averaged over 1000 realizations and keeping 
g = 1000 constant. It is worth noticing that the mean 
field prediction for a network of size N = 400 is that 
the system contains about eight communities on aver- 
age, which agrees well with the seven communities found 
by the local fitness approach used in Fig. [2j 

In addition we have calculated the average neighbor- 
hood degree as a function of node degree and found that 
our final networks show assortative mixing, i.e. high de- 
gree nodes are connected to other high degree nodes, as 
is common in other social networks. As a kind of hypo- 
thetical exercise we have found that if we modify Eq. ([5]) 
by substituting \xi— Xj \ with +2^ | and Eq. © by sub- 
stituting \xi + Xj\ with \xi — Xj\ thus favouring similar 
opinions in link deletion and opposite opinions in local 
rewiring, respectively, the final networks shows dissorta- 
tive mixing instead. 

Furthermore, though not reported in detail in this pa- 
per, we have also performed numerical simulations to in- 
vestigate the effect of y (i.e. the proportion of focal clo- 
sure events in Eq. @) on the final configuration of the 
system. We have performed calculations in networks of 
size iV = 400, keeping (ko) = 4 and g = 100 constant, 
and varying y from to 1. The main result is that the 
clustering coefficient remains practically constant, except 
when j/al, where it decays. The average shortest path 




FIG. 3: The ratio (k) /(N (C» as a function of g, for vari- 
ous values of N. Each point corresponds to an average over 
100 realizations. The inset shows the number of communities 
X as a function of N for g — 1000 and averaged over 1000 
realizations. 



is reduced when y increases, which is to be expected, 
since focal closure offers the opening of long jumps in the 
network in the same manner as disorder produces small 
world properties in a ring. Finally, the average degree 
increases slightly with y, also explainable by the same 
token. More studies in this direction will be made in the 
future. 



IV. MEAN FIELD APPROXIMATIONS 

In order to understand better the numerical results 
we analyze our model by using a mean field approach. 
With this one can investigate, for instance, how rapidly 
the agents reach an irrevocable opinion, that is, how the 
mean number of decided (or undecided) agents behaves 
as a function of time in early stages of the dynamics, 
or when t ss 0. If we take hi = as in our numerical 
calculations, Eq. ^ is symmetric with respect to the 
sign of xi. In this equation the term aifo({xj}) is likely 
to be very small for two reasons: 

1. for t « the term fo({xj}) ~ 0, because the distri- 
bution of Xj is symmetric around zero, and 

2. ai is a random parameter centered at 0, so when 
averaged over i it turns out to be very small. 



Then one could write 



iemi(i) 



which means that Xi decays or grows exponentially at 
a rate given by the magnitude of 7^ that varies locally. 
Without loss of generality we consider those Xi that go 
to +1. Of all the agents with Xi > only the ones with 



7 



7t > will reach eventually the maximum value 1, and 
this number is of the order of N/4. 

For early times, the number of agents with Xi = 1, 
denoted by N\ , is small and the probability of having a 
system of N\ agents in a "bath" of N should be of the 
canonical form Pi(Aq) = Nl /Z', where the "chemical 
potential" should fix the average number (Ni). Since the 
probability for all the states with a given value of X4 is 
of the same form, and N = Nq + N\, the probability of 
having N states with x ^ 1 is PoC^o) = er^ Na /Z, where 



z = E 



N o =0 



1 - e-v 



e 2 



sinh 



(^) 



(8) 



sinh (7^ 



Now, the average fraction of undecided agents (no) = 
(N ) /N at a given time can be calculated as 



1_ A No 
7, 2-, M ( 



No=0 



-fiN 



1 d\nZ 
"N dfi 



(9) 



The quantity n is difficult to calculate even with mean 
field approach so we will consider it as an adjustable pa- 
rameter. Here we assume that /x is of the form t/r, where 
r is a large number. Taking into account both values of 
sign(a;°) one finally gets 



N + 1 /N + lt 

(no) = 1 coth 

V 07 N \ 2 t 



1 , f t 
— coth — 

N \2t 



1 - B 



N 



(10) 

where P/v (t/r) is the Brillouin function found also for 
magnetic systems. In Fig. |4]we show the comparison of 
the mean field prediction of Eq. [TO] with numerical calcu- 
lations performed in networks of different sizes, and with 
a fixed value of g = 10 5 . The parameter r was fitted 
with a least squares technique and has a value of approx- 
imately gN/ TO, which for N > 50 is indistinguishable 
from the Langevin classical limit, independently of N. 
Then the characteristic time taken by an agent to reach 
its definite opinion is g/20, so almost all the agents have 
reached their definite opinions before rewiring. 

Notice that this result does not depend on the size of 
the network, and that the fitting is very good for short 
times, but starts to deviate significantly and systemati- 
cally for longer times. In the numerical calculations more 
agents remain undecided than what is predicted by our 
mean field approximation. This is to be expected, since 
we have not taken into account the effects of frustration 
due to changes of opinion produced by the distribution 
of the personal attitude parameter a. Indeed, the differ- 
ence between the mean field prediction and the numer- 
ical value can be used as a measure of the amount of 
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FIG. 4: Average fraction of undecided agents (no) as a func- 
tion of time, in units of g, for various networks of different 
size, where each point is an average over 1000 realizations. 
The mean field prediction fitted with Eq. (|10[) is shown in 
black. 



frustration, which can become important in subsequent 
rewirings. 

So far we have analyzed the fast dynamics of the state 
variable but have not considered the rewiring processes. 
An interesting thing to do is to compare the deviations of 
our model results from the known properties of a random 
network, for instance 



(k) = (C)N. 



(11) 



Notice that in our model (k) should be constant to first 
order, since the network changes are such that every time 
a link is deleted another one must be created. Therefore, 
a change of degree Afcj is due to simultaneous actions of 
both agents connected by a link in the parallel rewiring 
scheme. 

In each rewiring event every agent deletes and creates 
the same number of links, so the agent's degree does not 
change. However, if a neighbor decides to delete the same 
link and create a different one, then two links are created 
but only one deleted resulting in the average degree of 
the network to increase. Similarly, if two agents inde- 
pendently decide to create a link between them and both 
have deleted different links, the average degree decreases. 
This can be expressed as follows 



Ak, = Afc„ 4 



Afc; 



(12) 



where the first term stands for the correlated processes 
associated to link deletion increasing the degree and the 
second term for the processes associated to link creation 
decreasing the degree. Of the fcj first neighbors of node 
i that could be deleted, only a fraction / 2 + of them will 
produce a correlated increase of the average degree. Then 
one could write Afc,^ = fcj. 

In our numerical simulations we have considered only 
triadic closure mechanism as the process of creating 
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links. Thus, the second term in Eq. (|T2|) should be pro- 

(2) 

portional to the number of second neighbors n\ \ or 

Ak~ — f~n\ ' . Now, only a fraction of the (kj — 1) 
second neighbors of i connected to first neighbor j will 
produce a correlated decrease of the average degree. Let 
us denote this number by {fjkj — 1). However, this 
holds only in a tree, where all the second neighbors are 
different, but in a general network there should be cyclic 
closures, like triangles (n A ) and squares {n a ). Then one 
may write 



,( 2 ) 



E 

jGmi(i) 



(/r*j 



l)-2n 



A 



(13) 



where the summation runs over the first neighbors of i. 
The factor of 2 in the second term is due to the fact that 
there are two triangles associated with a single triad. The 
exact count of links, according to Eq. (TT2"]) is 



(Ak) = 



1 



N 

N E 

i=i 



fth-Sl 



E (/; 



k 3 - 1) - 2< 



U'emi(i) 



(14) 



The number of triangles is related to the clustering coef- 
ficient as 



C, 



2n 



ki(ki 1) 



(15) 



So far we have not made any approximations. The 
factors /j + and /~ are time dependent through the dy- 
namics of the state variables X{. It is clear that these 
factors are non-zero only if the two sites involved in the 
same rewiring process are still undecided, that is, the 
probability of having one of these correlated rewirings 
should be proportional to (no) , and can be calculated 
numerically by keeping track of the rewiring matrices Py 
and Qij defined as 



P, 



othewise 



Pij > 



and 



if jeivf !^ >0 

othewise 



Qij — 



where and 



(16) 



(17) 



are subsets of the first and second 
neighbors of agent i, such that their cardinalities are the 
same because rewiring is the only operation allowed (see 
Eqs. © and ©). Then it is clear that 



f + = 



(P 2 )i 



(P 2 h + {Q 2 )t 



and /. = 



{Q 2 )i 



(P 2 h + (Q 2 )* 



(18) 



The normalization factor comes from the fact that only 
operations taking place by both agents i and j simul- 
taneously contribute to the count in i. This is pro- 
portional to 5Zj(Py + Qij) 2 where PijQji = 0. Then 
/r = 1 — /+ = fa. One can neglect the term n? in 
Eq. pT)l because the number of squares should be small 
as compared with the number of triangles. Also, we may 
assume that fi is very similar in all sites and near to 
1/2. In fact one can trace the mean value of fi as a func- 
tion of time in the numerical calculations and see that 
it converges very rapidly to the value 1/2. Performing 
the average over i and assuming the degree, the fraction 
of nodes, and the clustering coefficient statistically inde- 
pendent we obtain 

Afc - (n ) 2 [(1 - })k - fk(fk - 1) + fCk(k - 1)] , 

where we have suppressed the brackets indicating av- 
erages and dropped the subindex in all the quantities. 
Therefore, we find a fixed point when 



k = 



1-fC 

f{f-cy 



(20) 



In Fig. Owe plot the result of Eq. (gDJ) with / = 1/2 
(continuous line) and compare it with the numerical re- 
sults (open circles) and the linear dependence predicted 
by Eq. (fTTj) for a random network (dashed line). Note 
that for any given average degree k the clustering coef- 
ficient for the mean field and numerical results is larger 
than for a random network, and for k — N/2 all the 
clustering coefficients meet at value 1/2. Also we would 
like to point out that although in the mean field approx- 
imation we have taken only a part of the second order 
effects into account, it behaves qualitatively the same as 
the numerical result. 



400 



350 ■ 
300 • 
250 ■ 
200 ■ 
150 
100 
50 



O N = 400 
• — random 




FIG. 5: The average final degree as a function of the average 
final clustering coefficient obtained from simulations in net- 
works of size TV = 400 and g — 100, and several initial degrees 
(ho) (open circles). The prediction by Eq. for a random 
network is shown as a dashed line, and the mean field result 
of Eq. (|20p is shown as a continuous line. 
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V. CONCLUDING REMARKS 

In this paper we present a general approach to the 
problem of opinion formation on an adaptive network, 
which can be considered as a basic problem of social 
network formation and communication. The dynamical 
equation Eq. (Q]) and its modification Eq. ^ contain sev- 
eral aspects of this problem, which to our knowledge have 
not been considered or addressed before. First, the role 
of general opinion of all the agents in the system repre- 
senting the overall mood sensed by each individual with 
his/her own attitude towards it. Second, the effect of the 
media as a kind of external field term not as a constant 
field as suggested before (see, e.g., [20] ) but as an agent 
dependent reaction term corresponding to a personal at- 
titude towards the "manipulation" attempts by the me- 
dia. This makes the system reminiscent to a random field 
model. Third, in our model the two time scales of the 
basic dynamical processes are clearly separated and de- 
tcrministically controlled, unlike in other studies where 
an additional parameter like bounded confidence controls 
the range within which different opinions act on each 
other [3Cj , thus affecting the time scale at which opinions 
change. We note that also in our model the bounded con- 
fidence parameter could be included, but we chose not to 
do so in the interest of not increasing the number of ad- 
justable parameters in the model. Fourth, in our model 
the rewiring has been done deterministically favoring the 
triadic closure mechanism rather than random rewiring 
mechanisms used in other studies [H, [H, [3(| [35| • In our 
model both the basic link-rewiring mechanisms of net- 
work sociology, namely focal and triadic closure, have 
been included but in this paper we have mainly focused 
on the effects of the latter rewiring mechanism. 

In this study rather than exploring the entire richness 
of our model, we have concentrated on the effect of the 
separation of time scales for comparison with previous 
results. We have shown that the important feature of 
coevolution is the separation of the two basic time scales, 
namely, the rapid dynamics of the state variable, and 
the slow dynamics of the network rewiring. The only 
parameter of our model has been <?, serving as a measure 
of this separation. 

As one of the key results we have found that for inter- 
mediate values of the time scale variable g the network 



turns out to organize in well-connected small commu- 
nities of agents with the same opinion. This is in accor- 
dance with the earlier study [3(| , where it was found that 
communities could form when the tolerance or bounded 
confidence parameter was varied, which could be under- 
stood to be due this parameter changing the effective 
time scale for rewiring. In addition we have studied fur- 
ther the role of the difference in time scales for fast and 
slow dynamics analytically by devising a mean field treat- 
ment. Using random network as a reference, we have 
found that the mean field results compare quite well with 
the numerical results and both of them differ significantly 
from the random network results. 

In the future we plan to investigate our model sys- 
tematically and in detail for the effects of personal atti- 
tude parameter a, and the personal random field term 
of an agent hi, as well as the relative importance pa- 
rameter y between local rewiring with triadic closure 
and global rewiring with focal closure mechanisms. Fur- 
thermore, another interesting point to investigate in our 
model would be the relative importance of short and long 
range interactions. 

As a final remark we believe that our fully dynami- 
cal approach and the kind of coevolving network model 
could be applied to various other situations, not only in 
social networks but in other fields, like symbiotic rela- 
tions between two species sharing an ecological environ- 
ment. Some of these applications are currently a matter 
of further explorations. 
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